Does the attractive Hubbard model support larger persistent 
currents than the repulsive one ? 
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Abstract 



We consider a one-dimensional Hubbard model in the presence of disorder. 
We compute the charge stiffness for a mesoscopic ring, as a function of the 
size L, which is a measure of the permanent currents. We find that for finite 
\q • disorder the permanent currents of the system with repulsive interactions are 

larger than those of the system with attractive interactions. This counter 
intuitive result is due to the fact that local density fluctuations are reduced 
in the presence of repulsive interactions. 

+- » 

a 

a 
o 

o 
> 

•i-H 



Typeset using REVTeX 



I. INTRODUCTION 

A hallmark of mesoscopic systems is the presence of permanent current upon applica- 
tion of an external fluxfl. Although a noninteracting electron theory is quite successful in 
describing qualitatively the features of such an effect, it strongly underestimates the magni- 
tude of the observed permanent currents compared to the observed experimental valuesuij. 
One possible way out of this discrepancy would be to include the electron-electron interac- 
tions. In general the combined study of disorder and interactions proves difficult to tackle, 
so that one has to resort to various approximations. In three dimensions Hartree-Fock-like 
calculations^ tend to suggest that permanent currents are enhanced by the interactions. Ar- 
guments based on level statistics reach the same conclusionsQu. In one dimension, on the 
other hand, numerical and analytical techniques alike allow to treat the interactions exactly 
enough so that firm results can be obtained. Surprisingly, however, various independent 
calculations performed on a spinless Hubbard model with nearest neighbor interactions3H2l 
or long range coulomb interactional have reached the opposite conclusion, namely that re- 
pulsive interactions are detrimental to permanent currents, which, therefore, casts doubts 
on the validity of the approximate calculations in higher dimensions. 

We here consider the full problem of interacting electrons with spin, in the presence of 
disorder, and examine the permanent currents in such a system. We will mainly concen- 
trate for the sake of simplicity on a purely local interaction (Hubbard model), but will also 
consider briefly finite range interactions. The interactions are treated exactly and we use 
a renormalization group calculationEjO to take care of the disorder. We show analytically 
that the permanent currents are enhanced by the presence of repulsive interactions, and 
that the result of the spinless model, although correct, was an artifact. As a byproduct of 
the study we also give a very simple formulation of the Bethe-Ansatz equations for the pure 
attractive Hubbard model. 

The plan of the paper is as follows. Section |Tj relates the so called charge stiffness to 
the permanent currents, and discusses various peculiarities that can occur for finite temper- 



ature. Section |(J introduces the model and treats the interactions using the bosonization 
procedure. This section is merely to fix the notations. Section [TV] computes the stiffness 
using a renormalization group calculation. Section |V| focuses on the case of attractive inter- 
actions, both with a large —U expansion and using Bethe ansatz, and section [VI] examines 
an extended (nearest neighbors) Hubbard model. Finally the conclusions of our study can 
be found in section |VT 



II. STIFFNESS VERSUS MEISSNER 

Instead of computing the permanent currents it is simpler to focus on the so-called 
conductivity stiffness -DE-i, which measures the strength of the Drude peak in a macroscopic 
system a{uj) = DS(uj) + cr rcg . The stiffness D can be related to the change of the energy of 
the ground state of the system in presence of an external flux by 
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Eq being the ground state energy of a ring in a field. Here <fi denotes the boundary angle 
= 27r///o where / is the flux threading the ring and /o = hc/e is the flux quantum. On 
the other hand, for a mesoscopic system, the permanent current measures the response to a 
finite flux by 



(2) 

Therefore the stiffness D provides a measure of the permanent currents for small (or close to 
a multiple of 2tt) flux since J = 2D(f). Although the complete calculation of the permanent 
currents at finite flux is also possible for a one dimensional interacting system, the calculation 
is more complicated in the presence of disorder, and the stiffness carries enough information 
for our present purposes. 

In order to compute the permanent currents, one should in principle distinguish between 
an odd and even number of electrons, assumed spinless for the present argument. For an 
odd number of electrons the energy is minimum in zero external flux, and the number of 
right and left moving electrons is equal. For an even number of electrons, due to the k — 
state, the number of right and left moving electrons cannot be equal, and the energy is not 
minimum at zero flux. If one has non-disordered noninteracting electrons, it is obvious that 
in that case the minimum of the energy corresponds to half a flux quantum. An external 
flux can be incorporated in the wave function by making the usual gauge transformation 

Ci -> e iAr % (3) 

where A is the vector potential. This transformation modifies the kinetic energy term and 
absorbs an external flux at the cost of a twist in boundary conditions. Since this transfor- 
mation is purely local and the interactions and disorder do depend on density only, it will 
not affect the interaction and disorder terms. Therefore the permanent currents for a system 
with an even number of electrons, even in the presence of disorder and interactions, will be 
identical to those of a system with an odd number of electrons, with a shift corresponding 
to half a flux quantum. In the calculation of the stiffness one implicitly assumes that the 
equilibrium state is a minimum of energy with zero current (equal number of right and left 
carriers). Therefore the stiffness for the case of an even number of electrons, measures the 
permanent current produced as a function of (0 — <po/2). 

The above connection between the permanent currents and the stiffness is valid only at 
zero temperature. At finite temperature the permanent current is given by the derivative of 
the free energy. One could think that the stiffness would be given by the second derivative of 
the free energy, within the canonical ensemble. But such a quantity is the so called Meissner 
fraction 

Ld 2 F 
Ps = 2W (4) 

Although both of p s and D are related to the current correlation function, they correspond 
to different limits. p s , being a thermodynamic quantity corresponds to the limit uo — > first 
and then to q — ► 0, whereas D which is related to transport corresponds to the limits taken 



in the reverse order. For finite T the two quantities are distinct. For a macroscopic system 
p s measures the superfTuid density and will be zero for a non-superfluid system, whereas D 
can be non zero if the system is perfectly conducting but not superconducting (e.g. free 
electrons in the absence of impurities). If one has a finite system then p s needs not to be 
zero even if the system is not superconducting. In general the Meissner fraction is given by 

L d 2 F _. .d 2 H. , 0Ei ^ e -PE a _ e -p E d H. , |2 .. 
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where n denotes a state of the system, H is the Hamiltonian and E n the energy of the 
eigenstates. We have separated out the contribution from the manifold of degenerate states, 
which are of especial importance for charge transport. In (|5|) the limits L — > oo and T — > 
do not commute in general. If one takes the limit L —>■ oo first, then p s tends to zero unless 
the system is really superconducting, as can be checked by computing explicitly p s for e.g. 
free electrons. For a mesoscopic system, where L is finite, p s will be finite and gives the 
slope of the permanent currents with the flux J ~ p s 0. One now has to distinguish whether 
we have an odd or even number of electrons, or more generally whether the ground state 
is not degenerate (odd number of electrons) or whether it has a twofold degeneracy (even 
number of electrons). In the first case, for small enough temperatures, the sums in fl5|) will 
be dominated by the ground state, other terms being suppressed by factors like e _/3A , where 
A would be the gap between the ground state and the first excited states. Such a gap 
remains finite for a systems of finite size. If the ground state is not degenerate then the 
average value of the current in zero external flux is zero in the ground state. One has 

pjTT 

(0|^|0) = (0| J|0) = O (6) 

The last term in (|5|) vanishes and one finds a positive p s , roughly temperature independent 
(dominated by the first term in (|5|)). There is, therefore, a paramagnetic permanent current 
for small flux. 

On the other hand, if the ground state has a twofold degeneracy, which occurs e.g. for 
the case of an even number of electrons, then in each of the degenerate ground states 4> v the 
current can (and will in general) be nonzero 

(0,1^10,) = (0,| J|<M^o (7) 

In that case (||) is dominated by the last term, which gives rise to a Curie type behavior, p s ~ 
— 1/T. There will therefore be a diamagnetic permanent current, with a slope proportional 
to 1/T. 

Many of these properties are well known for free electrons (see e.g.H3 and references 
therein) but the arguments presented here show that they are much more general and rest 
only on the degeneracy of the ground state and are valid for interacting electrons as well. 



III. ONE DIMENSIONAL HUBBARD MODEL: NOTATIONS 

Only a short derivation will be given here in order to fix the notations. More details can 
be found iiJH0. 

Let us consider for example the discrete ID Hubbard model with L sites 

h = -t Y, c l<r c j,v + U J2 n ^ n u ( 8 ) 

where () stands for nearest neighbors. Using the well known boson representation of fermion 
operators in one dimension the complete Hamiltonian (Q) becomes, away from half-fillingEjO 

H = H p + H a + -^-Jdxcos(V8(f) a (x)) (9) 



where H p and H a are defined by 

H u = — / dx 

2vri 






(10) 



IT and cf) are canonically conjugate variables and a is a short distance cutoff that can be 
identified with the lattice spacing. The p and a parts of the Hamiltonian (||) respectively 
describe the charge and spin degrees of freedom of the system. The gi± term is the scattering 
between electrons of opposite spins with an exchange of momentum of 2kp. The Hamiltonian 
([|) describes in fact the most general one dimensional Hamiltonian with spin conserving 
interactions, provided that the proper values of the K and u parameters are used. K p controls 
the charge excitations. K p > 1 corresponds to dominant superconducting fluctuations, 
whereas K p < 1 corresponds to charge- or spin-density wave (depending on the spin part of 
the Hamiltonian) dominant fluctuations. 

For the Hubbard Hamiltonian (Q) the various coefficients in (^) and ( |TUD are given, if U 
is small compared to Vp, by 

u p K p = u a K a =v F , u p /K p = v F + U/n (11) 

u a /K a = v F -U/n , gi± = U (12) 

For a macroscopic system gx± renormalizes to zero and K a to one for repulsive interactions. 
In the asymptotic limit u p , u a , K p are the only parameters needed to describe the long range 
properties and have been computed exactly using Bethe-AnsatzEa. For attractive interactions 
there is a gap in the spin sector, and only u p and K p are needed to describe the low energy 
properties of the model. They can also be computed from Bethe Ansatzt3 as will be seen in 
more details in section 0. 

The charge stiffness D can be obtainedli3c2l from the bosonized Hamiltonian ([]), and is 
simply given by D = 2u p K p . 

IV. EFFECT OF DISORDER ON D{L) 
Disorder can be added to (OF) by 



H dis = J dxW(x)p(x) (13) 

where W is a random potential. As is well known the forward scattering due to the potential 
does not affect the conductivity in one dimension (see e.g. [13]) and one can retain only the 2kp 
Fourier components of the random potential. Higher Fourier components are less effective 
since they do not scatter electrons on the Fermi surface, and do not correspond to low 
energy processes. A notable exception are Ak F components that will be discussed later. 



When expressed in term of the bosons variables (13) becomes 



H dis = j dx((x)e iV2 ~ Mx) cos(v / 20 (7 (a;)) + h.c. (14) 

where £(x) is a complex random potential corresponding to the part of the random potential 
having Fourier components close to 2k p. For simplicity we will take it Gaussian 

(ax)C(x')) = W c 5(x-x') (15) 

The effect of ( |T4"D can be computed using a renormalization group calculationt3't3, where 
one changes the cutoff (lattice spacing) a into ae l . We will just quote the results here 



dK p (l) 

dl 
dKjJ) 

dl 
du p {l) 

dl 
du a {l) 

dl " ~2 
dy{l) 




(16) 
(17) 
(18) 

-A(i) (19) 



dl 

dA(l) 



[2 - 2K a (l)]y(l) - A(l) (20) 

[3 - K p (l) - K a {l) - y(l)] A(/) (21) 



dl 

with the dimensionless quantities: 

The renormalizations of K v and u v are of first order in A and y 2 , and consequently can be 
neglected on the right sides of the first three equations. Contrary to the pure case, charge 
and spin degrees of freedom are now no more decoupled. 

For a macroscopic system, the physics implied by the equations (|T6| - |2T1) has been studied 
at lengthliBGJ. As can be seen from (|2"TD, there are two regimes depending on the initial 
values of K p , K a , y. A can scale to zero, the system is in that case delocalized, and as shown 
i n y,y j s dominated by divergent superconducting fluctuations. In the other regime A scales 
to large values, and the corresponding phase can be identified with a localized phase. In that 
case the RG equations break down below a certain length scale that can be identified with 
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the localization length. In the limit where A — > 0, as can be seen from (|2T|) the localized- 
delocalized transition occurs when 2 — K p = if y > (since y renormalises to zero and 
K a to one) or 3 — K p = if y < 0. For a mesoscopic system of size L, one can expect 
the size to play the role of an infra-red cutoff in the RG equations. When the renormalized 
cutoff ae l * ~ L, i.e. I* = \og(L/a), one can treat the disorder term in perturbation, and 
the stiffness is therefore simply given by the quadratic part of the Hamiltonian, with the 
renormalized coefficients D = 2u(l*)K(l*). Provided the size L is smaller than the length 
at which the equations cease to be valid one can use them to compute the various values of 
K, u, y as function of the size L of the system, and from that to obtain the stiffness D. Such 
a calculation is similar to the one performed for a macroscopic system to get the exponent 
K p c3 at finite temperatures. In that case the cutoff is provided by the dephasing length 
v F /T. 

The full dependence of the stiffness on the size of the sample needs a numerical integration 
of the RG equations, but the qualitative features can be understood by looking at the very 
small disorder limit. In that case one can neglect the renormalization of u,K,y in the 
equation for A, which gives 

A = A e (3_ ^" _A ' CT_9l±)log(L/a) (23) 

where A is the initial value of the disorder. Here we focus on the case of the Hubbard 
model. For small U one can use the values of the parameters ( O ) and one gets 

A = A (L/a) 1 - f//( ^ ) (24) 

Therefore the disorder grows more slowly for the repulsive Hubbard model than in the 
attractive one. In the same limit of a very small initial A, the stiffness is roughly given by 

D(l* = log(L/a)) = D(l = 0) - Cste / " A(Z) (25) 

Jl=o 

The bare stiffness D(l — 0) can be considered roughly independent of the interactions if the 
system is far enough from half filling as can be seen from ([TT| ) . The dependence in U comes 
only from lattice effectsBiiilScHl that breaks galilean invariance and are sensitive, for repulsive 
interactions, only for fillings close to a commensurate filling (mainly half filling where the 
pure system would be a Mott insulator). For attractive interactions the renormalization of 
the stiffness of the pure system due to interactions will become much more important and 
will be discussed in section [V} We will ignore in the following the change of the bare stiffness 
due to the interactions and will only focus on the effects due to the disorder. As can be seen 
from (|25|), the disorder term tends to drastically reduce the stiffness D, and this effect will 
be smaller for the repulsive model than the attractive one and the stiffness (the permanent 
currents) will be enhanced by repulsive interactions for a given size and a given disorder. 
For finite disorder one has to numerically integrate the RG equations. The result is 
shown in figure |I| and is in agreement with the simplified analysis above. This rather counter- 
intuitive result can be simply understood with the following argument: with the repulsive 
Hubbard model, the ground state is almost a spin density wave (with a power-law decay 
of the correlation functions) whose density is uniform. Such a ground state couples very 
weakly to non-magnetic impurities as is obvious from (|T3|). To couple to disorder, one has to 



distort the spin density wave and make a fluctuation of the density, a process that will cost 
an energy increasing with U. The disorder effect is therefore very weak, at least if the size 
of the system is not too big. On the other hand the attractive Hubbard model has a ground 
state that contains charge density wave fluctuations (although superconductive fluctuations 
are the dominant ones) which can get very easily pinned by impurities. On such a ground 
state the disorder will act very efficiently and drastically reduce the stiffness compared to 
the pure value therefore making the permanent currents smaller. Such an argument is in 
agreement with higher dimensions^. 

This is to be contrasted with a previously studied spinless modelBiiil. In that case both 
the attractive and repulsive ground state have density fluctuations, and both can be equally 
well pinned by disorder. Since in the attractive case the superconducting fluctuations tend to 
screen the disorder, the stiffness increases with attractive interactions. For the interactions 
to have a beneficent effect on the permanent currents one must necessary take a realistic 
model in which the main effect of the interactions will be to homogenize the density as is 
the case for the Hubbard model. 

If the interactions are infinitely repulsive, the system becomes equivalent to a model of 
spinless fermion with a Fermi momentum of 2kp. In that case, although the 2kp component 
of the disorder is inefficient (as is also obvious from the fact that it will no more correspond 
to a process on the new Fermi surface), one should worry about the 4kp component of the 
disorder. Such a Fourier component acts on the free spinless fermion, so that one recovers 
the stiffness of free electrons in the presence of disorder. The crossover between the two 
regimes would need a detailed analysis of the coupling of the 4/c^ component of the charge 
density to disorder which is way beyond the scope of this paper, but one could naively expect 
a maximum of the permanent currents for an intermediate value of the interactions. 

If the size of the system becomes large enough the disorder will renormalise to large values 
and the system will be localized. This is always the case for repulsive interactionsEjc3il3E3. 
For attractive interactions, a localized-delocalized transition is in principle possibleB^Ea't^lla 
(for K p > 3) and the stiffness could saturate to a finite value. We will show in the following 
section that for the particular case of the attractive Hubbard model, where one has only an 
on-site attraction, this transition does not occur and the system remains always localized. 

V. NEGATIVE U 

Let us consider the case of a large negative U. In that case one would naively imagine 
that the system should delocalize. In fact a very large on-site attraction cannot delocalize, 
and increases the localization (and therefore decreases the permanent currents). In the 
U — > — oo limit, one can perform a large \U\ expansion. Only pairs of particles can hop and 
if one introduces the operators bi = c^-j-c^, the attractive Hubbard model then becomes a 
model of hard core bosons with a hopping t' = t 2 /\U\ and a disorder A' = A. The residual 
interaction between the bosons is also on the scale of t 2 /U. In fact using the superexchange 
formulation or degenerate perturbation theory, the model maps on precisely to the 1-d 
Heisenberg antiferromagnet at a fixed magnetization (related to the density of particles), 
with an exchange energy 4t 2 /U. By a Jordan Wigner transformation this model corresponds 
to spinless fermion with a narrow bandwidth and with nearest neighbors interaction, in the 



presence of the old disorder. Such a system is obviously localized, and since the kinetic 
energy reduces with \U\ one expects the localization length to diminish when the attraction 
is increased. 

One can make the statements more quantitative for finite U by studying the RG equations 
in the attractive regime. In that case it is well known that there is a gap in the spin 
excitation spectrum and that only the charge sector remains ungapped. Keeping only the 
charge excitations into account the RG equations becomeli3't3 

dK (l) 1 9 . . , , . 

- -K 2 A(l) (26) 



dl 2 p 

du p {l) u 2 K f 



dl 

A(0 

dl 
With 



A(i) (27) 



[3 - K p (l)]A(l) (28) 



A(/) = (2C <y W e «)/W) ( 29 ) 

and C a is a constant of order unity coming from the <p a correlations in the perturbation 
expansion. The equations fllfiD can be used at scales above the size £ of a Cooper pair. This 
approach will therefore be adapted for reasonably large U. For small U it will be better to 
use the equations (|i~6-2~i~D, the crossover between the two regimes occurring when y ~ 1. 



In order to get the stiffness of the disordered system one needs the initial values of K p and 
Up in the absence of disorder as a function of the attraction U. As for the repulsive Hubbard 
model, they can be deduced from the Bethe-Ansatz solutionis. We will here give a derivation 
based on an appealing formulation introduced by Sutherland (in the absence of 0) involving 
the formation of Cooper pairs which scatter without diffraction. Let the number of particles 
be 2M, for which case we note the Bethe Equations for the attractive U Hubbard model, 
with energy E = — 4 J2j=i,m cos(Pj/2) cosh(Qj), where (Pj, Qj) = (Re, Im) arcsin(^ + i j), 
and ipj satisfy the Bethe EquationstZI 

2 
LP(ip j )/2 = 2irJj + 20 + J2 arctan[-(^ - ipj)] (30) 

i=l,M U 

where Jj are integers (half odd integers) for M odd (even). The flux comes in with a 
factor of 2 due to the charge of the Cooper pair. These equations can also be obtained 
from the repulsive casec3'E2l by using a particle hole transformation on the up electrons of 
a Half filled model, the spin excitations then map on to the above equations. This can be 
checked explicitly, using the idea of complementary solutions due to WoynarovichEj, which 
essentially rests on the recognition that the equations for real k's in the Bethe equations 
of the repulsive Hubbard model are the L real zeroes of a polynomial of degree L + IM. 
Hence the residue theorem of Cauchy helps in transforming equations involving the real k's 
to those over complex k's. The complex k's come in pairs, and are essentially pinned to be 
ipj ± ij, in order to satisfy the growth conditions. The error involved in writing down the 
above Cooper pair representation is of 0(exp(—L/£(U)), with £(U) the Cooper pair radius. 



The parameters u p and K p can be obtained by computing the compressibility x — u p /K p 
and the charge stiffness D = 2u p K p of the pure system from the Bethe ansatz ground state 
energy. Various values of u p and K p are plotted in figure |2] together with the stiffness D 
(for the pure system). As can be guessed from the large U expansion u p — > at large U. 
The fact that the parameter K p remains finite shows that the system remains interacting. 
One can check that the limiting value of K p is in agreement with the one obtained for the 
XXZ chained on which this system maps in the large U limit. Due to the reduction of 
the velocity, the stiffness of the clean system itself goes to zero at large U. Conversely 
to what happened for repulsive interactions where the stiffness of the pure system was 
nearly interactions independent, there is here a drastic reduction of the stiffness when the 
attraction is increased^. 

Using equations fllED one computes the stiffness in the presence of disorder. Here the 
main in the decrease of the velocity u p , which increases the relative strength of the disorder 
given by the dimensionless parameter fl2"5|). Some results are shown in figure ^, where we 
have normalized the stiffness to its value in the absence of disorder to avoid the trivial 
effect of renormalization of the bare stiffness by attractive interactions. In agreement with 
the previous section the reduction of the stiffness due to disorder becomes more and more 
important as the attraction U is increased. 

VI. EXTENDED HUBBARD MODEL 

In order to check the validity of the arguments presented here for a slightly more general 
model than the Hubbard model, we also look at an extended Hubbard model with a nearest 
neighbor interaction V defined by 

i 

In that case, and for small U and V the various parameters entering the equations are 

K p ~l — (2-cos(2£; F a)) (32) 

I'KVf KVf 

K <7 ~1 + + cos(2A; F a) (33) 

2txvf kvf 

U V 
y± ~ 1 cos(2k F a) (34) 

TTVf TTVf 

For small U and V the renormalization of velocities is of second order in U, V and can be 
neglected. When replaced in equations (PH) one gets 

(3-K p -K a -y) = l- ^— + — (1 - 2 cos(2A; F a)) (35) 

ZTlVp KVf 

If U 3> V the results are unchanged compared to the case of the pure Hubbard model. In 
order to check whether the physical ideas introduced here on the increase of the permanent 
current due to repulsive interactions are correct, or whether they are an artifact of the purely 
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local Hubbard model, one can consider the artificial limit where U = and V remains finite. 
Note that the model does not boil down in that case to the spinless fermion model since the 

V term still introduces interactions among opposite spins. 

In that case the effect of the interaction depends on the filling. For low filling, a repulsive 

V will tend to favor a spin density wave ground state again, whereas an attractive one would 
tend to pair particles on neighboring sites, giving a modulation of the density. For large 
fillings the situation changes: a positive V will now tend to favor two particles on the same 
site, to avoid paying the repulsion, and therefore to give a charge density wave. An attractive 

V favoring two particles on neighboring site will this time give a spin density wave. The 
change between a SDW towards CDW ground state occurs when K ff = 1 and as can be seen 
from (|33|) this will occur when kp = 7i"/4. Added to this is the competing effect that the 
more attractive V we have, the more there are superconducting fluctuations in the system 
which tend to reduce the disorder. Above quarter filling the two effects go hand in hand 
and repulsive interactions are detrimental to the stiffness, whereas below quarter filling the 
two effects will compete. From ( [35|) one can see that the point where a repulsive V again 
becomes favorable to the stiffness is kp = n/3. Below this filling, the fluctuations of the 
density generated by a repulsive V are too strong to be balanced by the superconducting 
fluctuations and a positive V will increase the permanent currents. 

VII. CONCLUSIONS 

We have looked in this paper at the stiffness constant of a Hubbard model as a function 
of the size of the system. The stiffness constant is directly related to the permanent currents 
in the presence of an external flux by J = D<p for small flux. We have shown that both 
the attractive and repulsive Hubbard model are always localized for a macroscopic system 
regardless of the strength of the interactions. In fact, using Bethe Ansatz solution or for 
large —U simple perturbation theory, one shows that the localization length decreases for 
attractive interactions due to the reduction of charge velocity. 

For a mesoscopic system, the stiffness in the repulsive Hubbard model is much less 
sensitive to disorder than for the attractive one. Therefore the permanent currents are 
enhanced by repulsive interactions. This surprising result is related to the fact that for the 
attractive Hubbard model the ground state contains strong charge density wave fluctuations 
that pins easily on the impurities, whereas repulsion favors a uniform density and makes 
the pinning harder. This property remains valid for a model with longer range interactions. 
In general the effects of the interactions on the permanent currents is controlled by two 
competing effects. One is the presence of density fluctuations in the ground state. The 
more there will be, the more easily the system will be pinned by disorder and the more 
the permanent currents will be reduced compared to the pure value. In general repulsive 
interactions will tend to favor a homogeneous density (local fluctuations in density will cost 
an energy increasing with the repulsion), and therefore will tend to increase the permanent 
currents. On the other hand, attractive interactions promote superconducting fluctuations 
in the system that tend to screen the disorder and therefore tend to increase the permanent 
currents. 

Previous studies of one dimensional systems, leading to the conclusion that repulsive 
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interactions reduced permanent currents, were performed on a spinless model. In such a 
rather artificial model, the first effect does not occur, since density fluctuations are always 
present both for attractive and repulsive interactions and therefore repulsive interactions are 
detrimental to permanent currents. In a more realistic model where the local interactions 
are the dominant ones (interactions in a real system do decrease with distance !), the density 
effect will dominate and the permanent current are increased. 

This study is, strictly speaking, restricted to one-dimensional systems, and a direct 
comparison of our results with experimental, three-dimensional, rings is not feasible. It 
nevertheless suggests that in mesoscopic systems the presence of repulsive interactions can 
considerably enhance the permanent currents, and confirms in the exactly solvable one- 
dimensional case, that the increase of the permanent currents is linked to a reduction of the 
local density fluctuations by the repulsive interactions. It is therefore tempting to ascribe 
the discrepancies observed between the measured and the computed (with a free electron 
theory) values of the permanent currents to such an interaction effect, an interpretation 
compatible with recent perturbative calculations^. 
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FIGURES 

FIG. 1. Normalized stiffness D/Dq as a function of the size of the system (in units of the lattice 
spacing a) obtained by numerically integrating the RG equations (|l~6|42~l|) . Dq is the stiffness in 
the absence of disorder. All energies are in units of the orignial Fermi velocity vf- The disorder 
W^/vf is fixed to W/vf = 5 1CT 4 . The full line is U/vp = 0, the dotted line U/vf = —0.5 and the 
dash-dotted line U/vf = 0.5. For a given size L, systems with repulsive interactions have a larger 
stiffness than those with attractive ones. 

FIG. 2. Values of u p and K p for the attractive Hubbard model as a function of the strength 
of the interaction \U\. These values are obtained by numerical integration of the Bethe-Ansatz 
equations for systems of L = 200 sites with respectively 90, 70, 50 particles per spin for the full, 
dotted and dash-dotted lines. This corresponds to density of n = 0.9, n = 0.7 and n = 0.5 particles 
per sites respectively. 

FIG. 3. Normalized stiffness D/Dq for the attractive Hubbard model as a function of the 
size of the system (in units of the lattice spacing a) obtained by numerically integrating the RG 
equations (26-ffiJ|). Dq is the stiffness in the absence of disorder. The effective disorder C a W^ 



is fixed to C a W^ = 5 10~ 4 and the density is n = 0.5 particles per site. The full, dotted and 
dash-dotted lines correspond respectively to \U\ = 5, \U\ = 10 and \U\ = 15. The corresponding 
bare stiffness are respectively Dq = 1.78, -Do = 1-04 and Dq = 0.72. Here again, for a given size L 
and fixed disorder the stiffness decreases with increasing attraction. 
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